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I.  INTRODUCTION 


One  of  the  dangers  of  using  electro-explosive  devices  (EED),  whether  for 
military,  commercial,  or  aerospace  use,  is  unintentional  detonation  of  the 
device  due  to  electromagnetic  energy.  In  order  to  evaluate  this  hazard,  several 
techniques  and  algorithms  have  been  developed,  mostly  for  single  frequency 
electromagnetic  sources. 

The  development  of  the  microcomputer  and  the  drastic  increase  in  the  Inci¬ 
dence  of  digital  implementation  which  followed  created  additional  problems  in 
the  area  of  electromagnetic  interference  (EMI).  The  interference  emitted  from  a 
computing  system  is  not  limited  to  a  single  frequency,  but  instead  falls  into  a 
broad  spectrum,  often  covering  several  hundred  megahertz  (MHz).  The  designer 
using  EED's  must  evaluate  the  power  contained  in  the  broadband  emissions  in 
order  to  insure  that  EED's  used  in  close  proximity  to  computing  equipment  do  not 
receive  enough  energy  to  accidentally  detonate.  This  report  presents  a  "quick 
and  dirty"  approach  to  evaluate  the  worst  case  energy  coupled  into  an  EED.  The 
actual  energy  received  will  always  be  less,  due  to  assumption  made  in  the  algorithm. 

The  equipment  required  to  use  this  method  consists  of  a  Hewlett-Packard 
8568,  a  spectrum  analyzer,  and  a  computing  controller  equipped  with  an  IEEE  488 
bus  interface. 

Communications  engineers  like  to  express  the  frequency  content  of  a  time 
domain  waveform  v(t)  in  terms  of  a  frequency  domain  spectral  density  g(f).  This 
spectral  density  is  effectively  a  measure  of  the  power  or  voltage  in  each  small 
frequency  component  of  the  signal.  The  total  power  contained  in  the  waveform 
may  be  found  by  integrating  the  spectral  p<wer  density  of  the  waveform  over  the 
frequency  range  of  interest.  Mathematically; 
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P  =  /  G(f )df  (1) 

fl 

If  the  waveform  frequency  domain  data  contains  a  series  of  discrete  fre¬ 
quency  components  of  amplitude  A,  then  integral  1  becomes  a  discrete  series  sum¬ 
mation. 
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If  the  bandwidth  of  the  receiver  is  chosen  carefully  so  that  the  power  contained 
in  one  small  increment  of  the  frequency  spectrum  does  not  "spill  over"  into 
other  increments  (Afs),  then  equation  2  allows  one  to  find  the  total  power  con¬ 
tained  in  each  increment  of  frequency. 

The  Hewlett-Packard  8568A  and  8566B  Spectrum  Analyzers  are  scanning  filter 
devices  which  presents  the  frequency  domain  data  as  a  1000  point  array  of  ver¬ 
tical  amplitudes  over  a  linear  range  of  frequency  steps,  so  that  Af ,  the  analy¬ 
zer's  frequency  step  satisfies  equation  3. 
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If  the  analyzer's  resolution  bandwidth  (RBW)  is  chosen  less  than  A f  (O.lAf 
or  less  is  ideal),  each  spectral  amplitude,  A  will  be  influenced  very  little  by 
the  signal  amplitudes  on  either  side  and  the  requirements  of  equation  2  are  met 
{f0  +  (M-l)  Af,  An_i  and  fQ  +  (M+l)Af,  An+j.}.  Since  the  HP8568A  and  8566B  allow 
the  user  to  preselect  the  RBW,  the  power  contained  in  a  broadband  waveform  may 
be  calculated  easily  by  summing  the  amplitudes  of  each  spectral  component, 
assuming  the  user  selected  the  RBW  properly.  This  may  also  be  accomplished 
under  software  control  of  an  IEEE  488  bus  compatible  calculator. 

The  system  this  lab  uses  is  composed  of  an  HP8528A  Spectrum  Analyzer 
connected  to  an  HP9825  Computing  Controller  using  the  HP  Interface  Bus  (HPIB). 

A  flowchart  of  the  program  used  is  shown  in  Figure  1. 

The  basic  algorithm  used  to  calculate  the  broadband  power  is  very  straight¬ 
forward. 

1.  Set  up  the  spectrum  analyzer  to  take  the  emissions  data  recorded  as  dByV. 

2.  Apply  a  frequency  dependent  correction  factor  for  the  antenna  effects  to 
convert  each  data  point  to  dByV/m 

3.  Convert  E  from  its  logrithmic  form  to  its  linear  form 
E(v/tn)  -  10 

4.  Convert  E  (v/m)  to  power  density  Pd(W/m2) 

Pd  -  E2/377(U/m2) 

5.  Calculate  the  worst  case  power  coupled  into  an  EED  (due  to  each  fre¬ 
quency  component)  assuming  no  losses  due  to  matching,  and  reception  of  signal 
due  to  a  resonant  dipole. 

p  -  Pd  (1«65)  X2  where  A  is  the  wavelength. 

4x 

6.  Sum  each  of  the  powers  due  to  the  individual  components  to  find  the 
total  power. 

II.  PROGRAM  DESCRIPTIONS 

This  program  is  divided  into  two  sections.  The  first  one  (lines  0  -  21)  is 
a  loop  which  controls  the  program  flow  by  calling  appropriate  subroutines.  This 
section  also  initializes  the  computational  and  data  acquisition  equipment.  The 
second  section  (lines  22-207)  contains  the  subroutines  called  by  the  control 
section. 

Each  of  the  subroutines  will  be  explained  in  detail  in  the  following  sec¬ 
tions  but  some  words  of  clarification  are  needed  concerning  the  control  section. 
Line  3  clears  the  HPIB  Interface  and  Line  19  stores  the  plotter  pen. 


Figure  1.  Flowchart 
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0:  *  15  NAY  81  POWER  SUMMATION  V 
1:  *  PQWSUN  •: 

2:  qsb  *D£V' 

3:  cli  ? 

4:  qsb  'PLINIT* 

5:  qsb  'ARRAY* 

6:  qsb  'qrat* 

7:  qsb  'TITLEl' 

8:  qsb  'title-iut* 

?:  asb  'SETUP* 

Id:  qsb  'FREQ* 

11:  qsb  'B8P' 

12:  qsb  'SWEEP' 

13:  qsb  'REDRL* 

14:  qsb  'REDSA' 

15:  gsb  'AHTCOR' 

16:  qsb  •SUHHATION* 

17:  qsb  'LEARN* 

18:  qsb  'STATEOUT' 

19:  oen* 

28:  dso  'END  OF  TEST '(beep; end 
21: 

22:  "ihhhhhhhhh»*«m**««  SUBROUTINES  «HHHHWW#*m*«**»***»***** : 
23: 

24:  "DEV  ASSIGNS  THE  HP1B  ADDRESSES': 

25:  "DEV‘: be?  •sa'!7t8!*lp*,6,'or'>785.‘25')719>*str*,714:r9t 


III.  SUBROUTINE  DESCRIPTIONS 
A.  DEV 

Subroutine  DEV  uses  the  HPL  device  statement  (DEV)  to  associate  string 
variables  with  peripheral  logical  unit  numbers  (HPIB  addresses).  Line  25  makes 
the  following  assignments: 


Device 

Name 

HPIB  Address 

HP8568  Spectrum  Analyzer 

sa 

718 

HP8966B  Line  Printer 

lp 

6 

HP9872A  Plotter 

gr 

705 

HP9825  Controller 

25 

719 

HP2631A  Line  Printer 

ptr 

714 

The  main  reason  for  using  the  device  statement  is  to  facilitate  program 
modification  if  the  HPIB  device  address  assignments  must  be  changed.  If 
peripheral  devices  were  referenced  directly  by  address,  changing  an  HPIB  address 
would  mean  modifying  every  statement  which  referenced  that  peripheral.  By  using 
the  device  statement  this  situation  is  averted,  and  only  the  device  statement 
making  the  assignment  must  be  changed  if  HPIB  addresses  are  changed. 
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B.  PLINIT 


This  routine  initializes  (clears)  the  9872A  Graphics  Plotter  and  sets  the 
plotting  window  scale  points  PI  and  P2.  It  also  sets  the  minimum  and  maximum 
scale  values  in  user  units. 

95:  •PLINIT*: 

96:  •INITIUZES  TIC  PLOTTER  AND  SETS  THE  VALUE  OF  PI  4  P2‘. 

97:  urt  WIN* 

98:  art  *gr*,*IP  1 88, 1808 ,13000  ,g008- 
99:  scl  -160,1100,-20,120 
181:  ret 

C.  ARRAY 

ARRAY  dimensions  the  necessary  storage  arrays  for  the  real  and  string 
variables.  The  HPL  dim  statement  is  equivalent  to  the  FORTRAN  DIMENSION  state¬ 
ment.  However,  there  is  one  difference.  In  HPL  the  string  variables  must  be 
dimensioned.  The  statement  dim  A$(2),  dimensions  string  variable,  2 
characters  long  to  be  named  A$. 

26:  ’ARRAY*: 

27:  ’ARRAY  DIMENSIONS  ALL  ARRAYS*; 

28:  tin  1UI88), 011881, LI1091).A*I2UM641,?$[B, 64] 

29:  ret 

D.  GRAT 

The  purpose  of  this  subroutine  is  to  draw  one  plot  on  which  the  final  data 
will  be  plotted.  Line  109  determines  the  size  and  tilt  of  the  labeling.  Lines 
111  through  118  draw  the  horizontal  lines  representing  the  linear  dBuV  scale. 

The  vertical  grid  lines  represent  the  linear  frequency  scale. 

107:  •grit*: 

108:  "draws  the  grand  on  the  plotter": 

109:  csu  1,2, 1,0 

118:  pen#  1 

111:  pit  8,8, -2 

112:  for  1*0  to  90  by  20 

113:  iolt  1088,0 

114:  iplt  0,10 

115;  iplt  -1090,0 

116:  iplt  8,10 

117:  neat  I 

118:  iplt  1800,0 

119:  pen 

120:  pit  0,0, -2 

121:  for  1*0  tt  400  by  101 

122:  iplt  0,100 

123:  iplt  100,0 

124:  iplt  0,-100 

125:  iplt  100,0 

126:  neit  I 

127:  iplt  I, 1*8 

128:  pen 

129:  ret 
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E.  TITLE1  and  TITLE-OUT 

TITLE1  inputs  an  array  of  string  variables  which  serves  as  a  graph  data 
identification  block.  The  input  format  is  eight  string  variables,  each  64 
characters  long.  TITLE-OUT  simply  prints  the  strings  input  in  TITLE1  on  the 
line  printer. 

84:  *TITL£1 

85:  ’INPUT  THE  DOCUMENTATION  INFORMATION  TO  BE  PLACED  ON  M  PLOTS’: 

86:  ent  ’STSTEM  NAME T',T$j ’SYSTEM  NAME;  ’4T$-*D$m ;  ”■»!$ 

87:  ent  ’TEST  DATE7 ■  fT* ; ’TEST  DATE:  ’6TVMM21;’’-tT$ 

38:  ‘MICOM/DRSNI-RTR’+Dtm 

89:  ent  ’YOUR  NAME? ’,Tt; ’NAME:  ’4T$-tDit4]:*’«T» 

90:  ent  ’TEST  NUMBER ?’,T$: ’TEST  NUMBER:  ’iT$-»0$[5!;’*->T* 

91:  ent  ’H0DE7\T4;’N0DE:  ’&THD«6];”-»T* 

92:  ent  •POLARIZATION?’,T$:’POLARIZATION:  *4TV>»[7)' ”->T* 

93:  ent  ’TEST  CONFIGUR AtlON? * , T» ;  "TEST  CONFIGURATION;  W-WfBl  ;”-»T* 

94:  ret 


130:  "Title-oet* : 

131:  ’ outputs  the  annotation  strings  entered  the  witor*' 
132:  art  *ptr’,char<10)4char<13) 

133:  art  *ptr *,char(27).i’r»’ 

134:  for  1*1  to  8 
135:  art  ’ptr’,D«I] 

136:  next  I 

137:  art  ’ptrVhar(27)i  V 
138:  ret 


F.  SETUP 


This  subroutine  sets  up  the  spectrum  analyzer  to  collect  the  data  using 
the  following  commands: 

IP  -  Run  the  spectrum  analyzer  instument  preset  routine.  Set  all  controls 
and  functions  to  their  default  states. 


KSC  -  Set  the  display  units  to  dBpV. 

ATODB  -  Set  the  spectrum  analyzer  input  attenuation  to  0  db. 


II  -  Select  Input  1.  The  HP8568A  has  two  signal  inputs.  Input  1  is 
selected  because  it  has  a  wider  bandwidth  (100  Hz  -  1.5  GHz)  than  input  2  (100 
KHz  -  1.5  GHz). 

30:  ’SETUP’: 

31:  ’SETS  THE  3568A  REFERENCE  LEVEL’: 

32:  art  ’so*, ’IP  KSC* 

33:  art  ’so’, ’ATODB’ 

34:  rtt 
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G.  FREQ 

This  subroutine  passes  the  start  and  stop  frequencies  for  the  band  to  be 
displayed  to  the  spectrum  analyzer. 


13i  -M8‘: 

34:  *FREQ  SELECTS  THE  FRE5UENCY  SPAN  CF  THE  ANALYZER*: 
37;  art  “sa‘, 'FB2S0HZ .FA20HZ* ;ret 


H.  BWP 

This  subroutine  passes  the  video  and  resolution  bandwidth  parameters  to 
the  spectrum  analyzer.  The  resolution  and  video  bandwidths  in  this  program 
represent  an  optimum  set  of  parameters  arrived  at  by  balancing  the  following 
requirements: 

I.  The  noise  floor  of  the  Instrumentation  system  is  required  to  be  at 
least  6  dB  below  the  MIL-STD  461  limits. 

2.  The  sweep  time  of  the  spectrum  analyzer  display  must  be  appropriate. 
Under  some  combinations  of  parameters,  the  sweep  times  may  be  as  long  as  25 
minutes;  much  too  long  to  run  repetitive  tests. 

38:  ‘BWP*: 

39:  ‘BWP  PASSES  THE  BANDWIDTH  PARAMETERS  TO  THE  3568A* • 

40:  w t  *Si', *RB10KZ,9BI0HZ* 

41:  ret 


I.  SWEEP 

This  subroutine  triggers  the  spectrum  analyzer  sweep,  using  the  A  trace 
memory  (trace  B  is  blanked).  The  trace  is  cleared  and  the  max  hold  (peak)  mode 
is  set.  The  spectrum  analyzer  is  allowed  to  sweep  continuously  for  approxima¬ 
tely  20  seconds.  This  helps  to  pick  up  transient  noise  with  a  low  duty  cycle. 


42:  'SWEEP*; 

43:  'SWEEP  CAUSES  THE  S568A  TO  SWEEP  THE  BAND  4  DISPLAY  IT*: 
44;  wt  *siVHD  Al* 

45:  wrt  *si*,*A2* 

46:  wt  *S4*,*S2* 

47:  wrt  ’si*(*TS* 

48:  wrt  *si*,*A3* 

49:  ret 
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J.  REDRL 


The  sole  purpose  of  this  subroutine  is  to  read  the  current  reference  level 
from  the  spectrum  analyzer  and  store  It  In  variable  R  for  later  use. 


1 01:  “REDRL  * : 

102:  "READ  THE  S563A  REFERENCE  LEVELS’: 
133:  w t  “si*  *RL  OA* 

104:  red  *54'  ,R 
105:  wrf  *54*,*Hf 
106:  ret 


K.  REDS A 

REDSA  is  the  subroutine  which  reads  the  raw  data  stored  in  the  spectrum 
analyzer  trace  memory  into  the  9825  controller  for  analysis. 

First  the  start  and  stop  frequencies  are  read  using  the  OA  (output  active 
function)  command.  The  start  frequency  is  stored  in  r2  and  the  stop  frequency 
is  stored  in  r3;  both  are  used  in  later  subroutines. 

The  trace  data  is  then  read  as  a  1000  element  array  of  vertical  points, 
with  the  element  number  representing  the  X  coordinate  and  the  vertical  magnitude 
representing  the  Y  coordinate.  The  displayed  magnitudes  range  from  a  minimum  of 
0  to  a  maximum  of  1040  in  screen  coordinates.  This  creates  a  square  field  of 
points,  each  side  of  which  is  1000  points  in  length,  which  leads  to  the 
Cartesian  coordinate  system  explained  above. 

The  FOR-NEXT  loop  contained  in  lines  61-64  is  the  portion  of  the  program 
which  transfers  the  raw  data  from  the  spectrum  analyzer  to  the  controller. 

Line  63  is  of  particular  interest.  The  first  portion  of  the  line  divides 
each  element  of  the  raw  data  array  (L[I])  by  1000  and  multiplies  by  10  to  con¬ 
vert  the  screen  magnitude  fo  dB  referenced  to  the  lowest  line  on  the  screen  gra- 
ticle.  The  next  portion  of  line  63  (R-100)  establishes  the  reference  level  for 
the  top  line  of  the  graticle  (using  subroutine  REDRL)  and  translates  that 
reference  level  to  the  bottom  line  of  the  screen  graticle.  The  two  numbers  are 
then  added,  producing  an  absolute  reading  in  dBuV. 

At  this  point,  the  data  stored  in  L[l]  should  correlate  with  the  raw  data 
read  on  any  EMI  receiver,  if  the  bandwidths  are  the  same,  the  data-point  fre¬ 
quencies  are  the  same,  and  the  instrument  tolerances  are  accounted  for. 


58:  “REDSA*: 

51:  “REDSA  READS  THE  TRACE  DATA  FROM  THE  8568A“: 
52:  “first  read  the  stir?  and  step  freogencies’: 
53:  fat 

54:  art  “si'/FA  OA* 

55:  red  *sa*,T:T-»r2 
56:  art  “sayFB  OA’ 

57:  red  ’si“,T;Hr3 
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58:  wrf  *sj 

59:  ‘now  read  the  data  in  screen  M.Y)  roordinstes*: 
68:  art  'sa'/Ot  TA* 

61:  far  1=1  t»  1001 

62:  red  "saVLlll 

62:  tUJ/10(IO*m+(IMOO)-»Lm 

t>4:  next  I 

65:  ret 


L.  ANTCOR  AND  ANTCORLIM 

These  subroutines  are  closely  related  and  will  be  considered  at  the  same 
time.  ANTCOR  takes  the  data  array  output  from  CHOP  or  BCOR  and  corrects  for  the 
effects  of  the  antenna  factor.  The  corrected  data  is  plotted  on  the  graph 
generated  in  GRAT.  Due  to  the  complexity  of  ANTCOR,  it  will  be  explained  in 
detail . 


Line  68  calculates  the  frequency  difference  between  sampling  points  by 
subtracting  the  start  frequency  (r2)  from  the  stop  frequency  (r3)  and  dividing 
by  the  number  of  points. 

Line  72  calculates  the  frequency  of  each  measurement  point  and  converts  to 
an  integral  number  of  MHz . 

Line  73  calls  ANTCORLIM,  which  supplies  the  slope  and  intercept  to 
generate  a  log-linear  function  used  to  approximate  the  frequency  currently  con¬ 
tained  in  rl.  (The  approximation  will  be  explained  in  detail  later). 

Line  74  Initializes  r6. 

Line  76  calculates  the  antenna  correction  factor. 

Line  77  adds  the  antenna  factor  to  the  data. 

Lines  79  and  80  actually  plot  the  corrected  data  [each  datapoint  of  L  (I)]  . 

Line  79  is  necessary  to  plot  the  first  point  in  each  range. 

Line  80  plots  all  others. 

ANTCORLIM  contains  a  set  of  variables  used  to  generate  a  linear  approxima¬ 
tion  to  the  antenna  factor.  These  variables  were  calculated  using  the  following 
procedure: 

66:  •wnCM': 

67:  ‘THIS  ROUTINE  CALCULATES  THE  FREQUENCY  DELTA*: 

68:  r3-r2/l®M-*M 

69:  pen*  l 

71:  *t»  slot  data* : 

71;  f«r  1=1  t«  HOI 

72:  <r2MI-l)«f4)/le6-*f 
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'3;  qsb  “ANTCORLIM* 

74:  3->r6 

”5:  ‘naw  qwerate  the  noianoiual  approximation  t,  .ne  3(,tenna  factor*- 
76:  3»l*qtF«le6MA-»r6 
77:  r6H.dl-H.il] 

~8:  if  1=1  ;plT  I,UIl.-2;ota  *2 
79:  if  1=1 pit  1,1111,-2  ' 

3D:  pit  I, LIU 
91 :  next  I 
82:  pen 
93:  ret 


1.  The  antenna  correction  factors  given  by  the  manufacturer  were  plotted 
on  semi-log  paper. 

2.  The  curve  on  the  log  chart  was  then  approximated  by  a  series  of 
straight  lines  over  frequency  intervals. 

3.  The  starting  and  stopping  frequencies  for  each  line  were  noted,  and 
the  slope  and  intercept  of  each  linear  segment  were  calculated  using  the 
following  relationships: 


S 


4  -r2 

log(r3)  -  log(rl) 


I  ■  r2  -  S  log  (rl) 


rl  ■  interval  starting  frequency 
r2  -  starting  amplitude 
r3  »  stop  frequency 
r4  =»  stop  amplitude 
S  *  slope 
I  *  Intercept 


These  slopes  and  intercepts  then  enable  the  calculation  of  the  antenna 
factor  for  any  datapoint  if  the  frequency  is  known.  The  frequencies  were  calcu¬ 
lated  in  line  72  using  the  3emilog  relationship: 


Y  ■  m  log  x  +  b 


139:  •ANTCORLIH*: 

140:  if  F>=2el  and  F<3el;-1.135774717-»S;l,929232528e!-'4:reT 

141:  if  F)=3*l  and  F<4*!;2.5625529eHS;-1.8*7081637e2-tA:ret 

142:  if  F>=4el  and  F<5d!;-2.W3770232eHS:1.708890511e2«A;reT 

143:  if  F)=5#l  and  F<6dl:-3.4J9898347el-*S;2,7452795«9p2->A;ret 

144:  if  F)=6el  and  F<7el;-1.94!84»b87el->S:1.6033?3056e2-‘A:reT 

145:  if  F>=7«1  and  F<8el;l .»34626412el-‘S;-7.316745634eHA;r»T 

146:  if  F>=8*1  and  F<9el;1.954937807e!-*S:-l,459tl04«41e2-'A:reT 

147:  if  F)*9el  and  F<le2;t  ,180134647e2-tS;-9,291577175e2-‘A:ret 

148:  if  F > = 1  e2  and  F< 1 . le2;-2. 657474371e! -*S ; 2 . 2?59,?9*99e2*A  r?t 
149:  if  F>=1  ,le2  and  FU.2t2;2,381673128eHS:-l  ,7761?6887e2-*A:r?t 
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150:  if  F)=1.2e2  and  F<1.3e2:.-2.876695653->S;3.30413A55?e!-tA;ret 

151;  if  F>=1 , 2e2  and  F<1.4e2;-S.835549619?HS;5.693326239p2-‘A'rpt 

152:  if  F)=l,4e2  and  F<1.5e2;6.807364332el-'S:-4.~68675?l?p2-’A;rpt 

153:  if  F>=1.5e2  and  F<1 ,6e2;4.638099184el-S;-3.o49l522!9e2-»A;r»t 

154:  if  F)=1.6e2  and  F<1.7e2;5.456773725el-*S;-5.141214634e2-»A'rpt 

155:  if  F>=t.7e2  and  F<1 . 3e2 : 2 . 81 9899046el -*S : -2 . ! 4790 3506e2^4 ret 

156:  if  F>=1 .8*2  and  F <  1 . 9e2 1-5 . 1 1 04940 19el -»5 ; 4 . 398852376 e2-»A ; r el 

157:  if  F>=1 ,9e2  and  F<2e2;9 . 875924846el -»S :  -8 . 0 C80  34838e2->A  ret 

158:  ret 

159:  'SUMMATION' : 

160;  0-»J 

161;  fet  6,38x,el2.5 

162:  fet  5 ,18i , **#*«»  W*»«SUNH  AT  !DN  OF  POWER (WATTS/'CQL'ARE  -ETER«*****«»" 

163:  for  1*1  to  1801 

164:  LUJ-tV 

165:  fntlV/20-6)->P 

166:  Pt2/377-*P 

167:  oto  +2 

168:  3e8/F-H. 

169:  P->D 
170:  gsb  ‘SUN* 

171:  next  I 
172:  urt  6,5 
173:  urt  6.6,J 
174:  ret 

175:  • *«******»SUBRCUT InES**»*******  * : 

176:  "SUM* : 

177:  8+J+J 
178:  ret 
179:  ■LEARN': 

180:  'THIS  SUBROUTINE  SAVES  THE  STATE  OF  THE  INSTRUMENT*; 

181:  die  C$111, 151 
182:  wrt  718,'OT* 

183:  for  N=1  to  11 
184:  red  718,C»tNl 
185:  next  N 
186:  ret 

187:  'STATEOUT': 

188:  for  1=1  to  11 
189:  prt  C$113 
190:  next  I 
191:  pit  -100, -5, -2 
192:  lbl  C$1103 
193:  pit  0,-10, -2 
194:  lbl  C$133 
195:  pit  450,-10,-2 
196:  lbl  C$14] 

197:  nit  725,-10,-2 
198:  lbl  C$151 
199:  pit  995, -5, -2 
200:  lbl  C$1111 
201:  pit  -100,95,-2 
202:  lbl  C$18] 
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233:  o It  -138,112,-2 
203:  lbl  ‘HP* 

285:  pit  -41,152,-2 
206:  lbl  C3I7J 
207:  Pit  240,102,-2 
208:  lbl  C»lb] 

209:  ret 
*6335 


M.  SUMMATION 


Subroutine  SUMMATION  obtains  the  power  contained  in  each  incremental  fre¬ 
quency  division  and  sums  them  over  the  band  of  Interest.  This  is  done  in  the 
following  manner. 

Line  160  sets  the  sum  counter  J  »  0. 

Line  161  is  a  format  statement  that  causes  the  line  printer  to  skip  38 
spaces  and  print  out  a  12  digit  number  with  5  digits  to  the  right  of  the  decimal 
point. 

Line  162  sets  up  the  format  for  the  label  “***  SUMMATION  OF  POWER  (W, 
SQUIBB  CORRECTED  ***". 

Line  163  and  164  interrogate  the  array  L(I)  where  the  1000  data  amplitude 
value  is  stored  in  the  variable  location  V(VOLTAGE). 

Each  successive  line  performs  the  following  algebraic  operations. 

Line  165  -  10  v/20-6»E 

Line  166  -  E2/377-P 

Lines  167-168  -  1.65  PL2/4it«D 

Line  169  -  turns  control  over  to  subroutine  ‘'sum”  which  adds  the  contents 
of  the  variable  D  contents  on  the  sum  variable,  I. 

Line  177  returns  control  to  line  170  which  recrements  the  loop  counter,  I. 

After  all  1000  power  amplitudes  are  summed,  lines  171-172  causes  the  line 
printer  to  write  out  the  power  summation  according  to  format  statements  5  and  6 
(lines  161-162). 

N.  LEARN 

The  subroutine  (lines  180-184)  dimensions  space  for  11  string  variables 
each  with  15  elements.  A  spectrum  analyzer  parameter  is  stored  in  each  string 
variable  location. 

O.  STATE0UT 

Lines  187-207  causes  the  plotter  to  print  contents  of  the  string  variable 
location  C$  on  the  output  plot. 
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